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Abstract-In this study, we develop new analytical and numerical 
speciation calculations corresponding to all kinds of 
combinations of complexation constants and hydrolysis constants 
given by more or less reliable literature data, and associated to 
fifteen probable uranyl-citrate species: M P L„, (with: M = (U0 2 ) 2+ , 
L = (HCit) 3 and p = 1, 2; n = 1, 2, 3), M p L„H m , (p = 1, 2; n = 1, 2; 
m = 1, 2, 4), MpCOH),, (p = 1- 4; q = 1-7), M p L n H m (OH) q , (p = 1, 2; 
n = 1, 2; m = 1, 2; q = 1, 2) as well as to the "improbable": 
[(U0 2 )H 2 (HCit) 2 ] 2 " specie (reported by only one reference), in 
order to analyse deeply their impact on the speciation diagram of 
U(VI) in citrate aqueous solutions at 25°C, in a wide range of pH 
and for citrate to uranium ratios R c = U(VI)/Citrate equal to 0.02, 
1 and 2. The main result of this work is that, although the 
employed combinations present some significant dispersions, 
they leads to the same result: the hydrolyzed monomer 
[(U0 2 )H(HCit)(OH)]-, the hydrolyzed dimer 

[(U0 2 ) 2 H 2 (HCit) 2 (OH) 2 ] 2 - and the monomer (U0 2 )(HCit)]" are 
largely predominant in the pH range of our interest (2< pH< 5.5), 
by comparison to the percentages of the other species 
presumably existing (except the [(U0 2 )H 2 (HCit) 2 ] 2 " specie). This 
investigation is therefore a supplementary proof of the reliability 
of the processes of protonation and hydrolysis of both the 
monomer and the dimer: [(U0 2 )(HCit)]", [(U0 2 ) 2 (HCit) 2 ] 2 " 
detected via the experimental electropherograms obtained by the 
capillary electrophoresis method which was used because of its 
capacity to separate the different hydrolysis products of the 
anionic U( VI) -citrate complexes. Consequently, it is thus possible 
to underline that the hydrolysis of U0 2 + which generally occurs 
at relatively elevated pH is facilitated at lower pH by the citrate 
complexation followed by the protonation. 

Keywords-Uranium,- Citrate,- Speciation ; Electrophoresis; 
Analytical Numerical 

I INTRODUCTION 

The ability of uranium to form strong stable complexes 
with organic compounds has been extensively utilized for 
cleaning up uranium contamination in soils and in the nuclear 
industry. Citric acid is a common constituent which forms 
strong water soluble complexes with U(VI) and U(IV) over a 
wide range of pH. The U(VI) citrate complexes have already 
been investigated in several experimental and modelling 
studies [1,2,5,6,10]. Accordingly, it was concluded that citrate 
has a strong affinity for uranium. However, it was also shown 
that this affinity varies with the experimental method, the 
ionic strength, the pH and the stoichiometry of the (M = U(VI) 
= U0 2 2+ )/(L = ligand = HCit 3 ") ratio. This dependence could 
partially explain the fact that the reported data by the literature 
concerning the different stability constants P pn for M p L n 
complexes scatter markedly. In fact, the M p L„ complexes 
could exist under different protonation and hydrolysis states: 
M p L n H m and M p L n H m (OH) q . However, most of the analytical 
techniques cannot differentiate between these different states. 



From recent analysis of all the previous works [8-10], a 
limited number of U (VI) citrate complexes have been 
selected ([(U0 2 )(HCit)]\ [(U0 2 ) 2 (HCit) 2 f, [(U0 2 )H(HCit)]) 
and their corresponding stability constants P pnm re-determined. 
More recently [16], we have underlined the advantage of the 
capillary electrophoresis technique over most of the speciation 
methods previously used because it should give the possibility 
to specify all the anionic M p L n H m (OH) q entities, by separating 
them on the basis of the differences in their electrophoretic 
mobility even at low concentration closer to the level found in 
nature. Accordingly, we have assigned some three peaks 
found in the absorbance spectra of our capillary 
electrophoresis experiments for pH < 5, to the following three 
species: the U(VI)-citrate monomer, the hydrolysis product of 
this monomer complex: [(U0 2 )H(HCit)(OH)]" and the dimer 
of the hydrolysis product: [(U0 2 ) 2 H 2 (HCit) 2 (OH) 2 ] 2 \ The 
latter two species were postulated in reference [2] and 
therefore play a crucial role in our speciation calculations 
despite the controversy concerning them from a chemical 
point of view. In particular, we have signalled for pH >3 the 
occurrence of an inverted tendency accompanied by a strong 
modification of the calculated speciation diagrams when we 
take into account these hydrolysed complexes via their 
corresponding constants: p lm and p 2222 . Indeed, the use of 
these additional constants leads to a completely different 
U(VI) distribution so that 60-70% of uranium is present in the 
form of hydrolyzed monomer rather than [(U0 2 ) 2 (HCit) 2 ] 2 ~ 
and 20-30% in the form of hydrolyzed dimmer rather than 
[(U0 2 )(HCit)]" Notice that the speciation diagrams given in 
our previous work paper (in the absence and in the presence 
of the hydrolyzed monomer and dimmer) are calculated on the 
basis of two constants: pn and p 22 selected among some 
published sets of values currently "more or less" established. 
However, in spite of the detection of the three peaks above 
mentioned, this inverted tendency is not quantitatively 
completely in agreement with the electropherograms results 
and we have explained this difference in assuming that the 
protonation and hydrolysis effects are non-reproducible 
kinetic processes. 

Instead, we will expose in the present work, new 
speciation calculations corresponding to all kinds of 
combinations of complexation constants given by more or less 
reliable literature data, and associated to fifteen probable 
species: M p L n , M p L n H m , M p (OH) q , M p L n H m (OH) q , as well as 
to the "improbable": [(U0 2 )H 2 (HCit) 2 ] 2 " specie (reported by 
only one reference), in order to analyse deeply their impact on 
the speciation diagram of U(VI) in citrate aqueous solutions. 
In the first step of this study we will therefore investigate 
theoretically (analytically and numerically) the different 
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speciation diagrams obtained at 25°C, in a wide range of pH 
and for different citrate to uranium ratios, following the use 
partly or completely of the uranyl-citrate complexes and 
uranyl- hydroxides mentioned above, and in assuming that the 
different reactions of complexation, protonation and 
hydrolysis have reached their complete chemical equilibrium. 
In the second step, we will compare the similarities and the 
differences between them in the perspective of verifying the 
fact that: for the pH range of our interest (2< pH< 5.5), despite 
their dispersion, the use of the several sets of values (more or 
less reliable), of the different complexation constants given by 
the literature, confirms the large predominance of the 
hydrolyzed monomer: [(UO z )H(HCit)(OH)]" and of the 
hydrolyzed dimer: [(U02)2H2(HCit) 2 (OH) 2 ] 2 ~ by comparison 
to the percentages of the other species presumably existing 
(except the [(U0 2 )H 2 (HCit)2] 2 " specie). This investigation 
could be therefore considered as a supplementary proof of our 
previous interpretation of the experimental electropherograms 
obtained by the capillary electrophoresis method (in particular 
the reliability of the processes of protonation and hydrolysis 
of both the monomer and the dimer: [(U0 2 )(HCit)]~, 
[(U0 2 ) 2 (HCit) 2 ] 2 ~), and also as a test of the validity of the 
different literature data as well as of the assumption of the 
chemical equilibrium. 

II METHOD OF CALCULATION 

A. Definitions 

We will designate in the following definitions by [X;] the 
molar concentration of specie i and by y X ; the activity 
coefficient of X; in the Henry reference in water at 25 °C 
which is related to the ionic strength I according to the 
generalized Debye equation [y] : 

log(yxi) = -0.5115(Zi) 2 I 1/2 /(l + 1.5I 1/2 ) + S jeij (Cj) (1) 



I = 0.5Z j (Z j ) 2 C j 



(2) 



where Zj and Cj are respectively the valence and the molar 
concentration of the j species, and e (/ is a specific interaction 
coefficient between species i and j. 

1 ) Complexation 

We define the complexation constant P' pn and the apparent 
complexation constant P pn corresponding to the following 
complexation equilibrium respectively by: 

pM + nL <^> M p L„ : P pn = [M p L n ][L ]" n [M]" p (3a) 

P'pn =(YM P Ln)(YL)- n (YM)- P Ppn (3b) 

2) Protonation 

We define the protonation constant P' pnm and the apparent 
protonation constant P pnm corresponding to the following 
protonation equilibrium respectively by: 

M p L n + mH + o M p L„H m : P pnm = [M p LnH m ] [M p L n T 1 [H , T m 

(4a) 

(4b 



P'pnm- (YMpLnHmXYMpLn) '(Yh + ) m Pp 



3) Hydrolysis 

We define the hydrolysis constant P' pnmq and the apparent 
hydrolysis constant P pnmq corresponding to the following 
hydrolysis equilibrium respectively by: 

M p L n H m + qH 2 o M p L n H m (OH) q + qH + : p pnmq = 



[M p L n H m (OH) q ] [H^tMpUHJ- 1 (5a) 

P'pnmq= (YMpLnHm(OH)q)(YH + ) q (YMpLnHm) 'P'pnmq (5b) 

a) Peculiar Case (n = 0; m = 0) 

We define the hydrolysis constant KH' pq of the metallic 
ion M and the apparent hydrolysis constant KH pq . 
corresponding to the following hydrolysis equilibrium 
respectively by: 

pM + q' OH " o M p (OH) q . : KH pq . = [M p (OH) q .] [OH] " 
q '[M]- p (6a) 

KH' pq = (YM P( OH) q 0( YoH)- q '(YM)- p KH pq . (6b) 
4) Acidity 

We define the acidity constant K a i and the apparent 
acidity constant K a i corresponding to the following 
equilibrium respectively by: 

L + 1H + <^HL : K a 1= [H^HUT^hV (7a) 



K a 'i = (YHiL)(YL)' 1 (YH + )"'K a ' l 



(7b) 



III EXPRESSIONS OF THE SPECIES CONCENTRATIONS IN TERMS 
OF [M], [L] AND [H+] 

A. General Expressions 

The expressions of the different species concentrations in 
terms of [M] and [L] concentrations are given by: 

[M p L n ] = p pn [M] P [L ] n (8) 

[M p L n H m ] = p pnm p pn [M] p [L ] n [H + ] m (9) 

[M p L n H m (OH) q ] = p pnmq p pnm p pn [M] p [L ] n [H + ] m [H + ]" q (10) 

[M p (OH) q .] = (KH pq .) [OH] q '[M] p (11) 



[H 1 L] = (K a 1 )[L][H + ] 1 



(12) 



B. Equations of Conservation 

If we define C° M and C° L as, respectively, the total 
concentrations of M and L, thus: 

C° M = [M] + Z P p[M] p {Z n (P pn + Z m p pnm p pn [H + ] m + Z m Z q 

Ppnmq Ppnm Ppn [HTlHTXL]"} 

+ Z P p[M] p { E q . (KHpq.) (Kw) q '[H + ]- qi } (13) 

C° L = E p [M] P {E n n( P pn + ImPpnm Ppn [HT + Imlq Ppnmq 
Ppnm Ppn [HT[HV )[L f } +{ 1+ I, (K 1 ,) H*] 1 }[L] (14) 

KW = [H+][OH-] is the apparent ionization constant of 
water. 

Now, if we note: x = [M] and y = [L], we will obtain the 
following equations: 

C° M = x +E p px p LA^y" +X P px p E+ {1+ Ei (K a ,) H + ]' } [L] 

(15) 

C° L = I p x p I n nA np y n + Dy (16) 

^Vith . A n p — ( Pp n + SmPpnm Ppn [H ] SmSq Ppnmq Ppnm Ppn 

[H + ] m [H + ]- q ) (17) 



E p = [E q . (KHpq.) (K w ) q '[H + ]" q '] 
D= {1+1, (K a i) [H + ]'} 



(18) 
(19) 
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Note that the parameters A np , E p and D depend on the 
different apparent equilibrium constants and are function on 
the ionic strength I = 0.5S i (Z i ) 2 C 1 and pH. We proceed 
therefore to a recurrent self consistent method: As a first 
approximation we equal the apparent constants to their 
corresponding thermodynamic constants, we calculate 
therefore for a given pH, the parameters A np , E p and D 
according Eqs. (17-19) in order to obtain a system of two 
equations with two variables: x and y according Eqs. (15-16). 
The resolution of this system for a given total concentrations 
C°m and C°l allows us to calculate the different concentrations 
of all the species present in solution: [M p L„], [M p L n H m ], 
[M p L n H m (OH) q ],[M p (OH) q .] and [HL] according Eqs.(8-12) 
and therefore all the activity coefficients via the ionic strength 
I. In the second step we recalculate the new apparent 
constants via Eqs. (3b -7b) and we repeat the procedure until 
convergence. 

The resolution of the algebraic system at each step is 
achieved using a basic numerical program involving n 
successive buckles and verifying the convergence 
conditions: [(x„ - x„_i) /x„] <10" 5 and [(y n - y„-i)/x n ] <10" 5 . This 
method allows us to establish for a given C° M and C° L , the 
corresponding speciation diagram giving the variation of the 
fraction in % of the different species with the pH. 



IV THEORETICAL URANYL-CITRATE SPECIATION DIAGRAM IN 
AQUEOUS SOLUTIONS 

Now, we will apply this general method in the case of: 

M = (U0 2 ) 2+ , L = (HCit) 3 ", in order to obtain the speciation 
diagram of Uranyl-Citrate complexes in aqueous solutions. 
For this, we will use different literature data for the 
calculation of P pnmq , P pnm , P pn , (K p H q ) and (K a ,) constants. 

A. The Most Probable Uranyl-Citrate Complexes in Aqueous 
Solutions 

The most probable uranyl species present in citrate 
aqueous solutions, according to the literature [2,3,7,10], are 
collected in Table 1. 

B. Calculation of Ppnmq, fipnm and fipn Constants from 
Literature Data 

In order to calculate P pnmq and P pnm constants, we 
proceeded as follows: 

According to the recent and less recent literature data 
[2,3,7,10], we can calculate the four protonation constants 
(Pin), (Pi 12) (P122) , (P224) and the two generalized constants 
(Pmi) and (P2222), in terms of the intermediate constants (Kj), 
(K 2 i), (K 22 i), (K 22 4i), (Kn) and (K m ) corresponding to the 
equilibria given in Table 2. 



TABLE 1 THE DIFFERENT URANYL SPECIES FORMED IN CITRATE AQUEOUS SOLUTIONS 



Complexes: Ppn 


Protonated complexes: fS pnm 


Uranyl 
hydroxides: KH pq . 


General hydrolyzed complexes: p pnmq 


[(U0 2 )(HCit)r 


[(U0 2 )H(HCit)] 


[(U0 2 )(OH)] + 


[(U0 2 )H(HCit)(OH)] 


[(U0 2 )(HCit) 2 ] 4 - 


[(U0 2 )H 2 (HCit)] + 


[(U0 2 )(OH) 2 ] 


[(U0 2 ) 2 H 2 (HCit) 2 (OH) 2 f 


[(U0 2 ) 2 (HCit) 2 f 


[(U0 2 )H 2 (HCit) 2 f 


[(U0 2 )(OH) 3 r 




[(U0 2 ) 2 (HCit) 3 f 


[(U0 2 ) 2 H 2 (HCit) 2 ] 


[(U0 2 ) 2 (OH)] 3+ 






[(UQ 2 ) 2 H 4 (HCit) 2 ] 2+ 


[(U0 2 ) 2 (OH) 2 ] 2+ 








[(U0 2 ) 3 (OH) 4 ] 2+ 








[(U0 2 ) 3 (OH) 5 f 








[(U0 2 ) 3 (OH) 7 r 








[(U0 2 ) 4 (OH) 7 ] + 





TABLE 2 HYDROLYSIS AND COMPLEXATION CONSTANTS OF URANIUM AT 25° 



Equilibrium 


log(K) : Ref [2] 


log(K) : Ref [10] 


(U0 2 ) 2+ + H(HCit) 2 <=> [(U0 2 )H(HCit)] 


log(K,) = 6.0 (1 = 0.136) 


log(Ki) = 4.23 (I = 0.1) 


(U0 2 ) 2+ + H 2 (HCit) <=> [(U0 2 )H 2 (HCit)] + 




log(K 2I ) = 2.79 (1 = 0.1) 


(U0 2 ) 2+ + 2[H(HCit)] 2 o [(U0 2 )H 2 (HCit) 2 ] 2 " 




log(K 22I )= 11.2(1=?) 


2(U0 2 ) 2+ + 2[H 2 (HCit)]" o [(U0 2 ) 2 H 4 (HCit) 2 ] 2+ 




log(K 224 ,) = 8.9 (1 = 1) 


(U0 2 ) 2+ + [H(HCit)] 2 + H 2 
[(U0 2 )H(HCit)(OH)] + H + 


log(Kn) = 2.84 (1 = 0.136) 




2(U0 2 ) 2+ + 2[H(HCit)] 2 + 2H 2 o [(U0 2 ) 2 H 2 (HCit) 2 (OH) 2 ] 2 + 
2H + 


log(Km) =7.68 (1 = 0.136) 
and 9.04 (I = 0.) 
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The different expressions of the P pnm and P pnmq constants 
are summarized in "Table 3". 

TABLE 3 EXPRESSIONS OF BPNMQ AND BPNM CONSTANTS OF URANYL- 
CITRATE COMPLEXES 



Pill 


Pun 


Pi 12 


(K'NXK.MP,,)- 1 


(KnXK,) 1 


(K-jXKjOCP,,)- 1 


P122 


(3222^2222 


^224 


(K a i) 2 (K 22I )(Pi 2 )' 


(K III )(K a 1 ) 2 (P22)-' 


(K a 2 ) 2 (K 224 i) (p 22 )-' 



C. Equations of Conservation in the Case of the Uranium 
Complexation with Citric Acid 

According to the data given in Table 1 , the most probable 
values for p, n, m and q are: p = 1, 2; n = 1,2 ; m = 1,4; and q 
= 1,2 ; the corresponding parameters A np and E p are thus: 

A n = A = (pn+Piiipii [H + ] + P„ 2 p„ [H + ] 2 +P1111P111P11) 

(20) 

A 21 = (P 12 +p 122 p 12 [H + ] 2 ) (21) 
A, 2 = p 2i ( is neglected) (22) 
A 22 = C = (P 22 + p 222 p 22 [H 1 "] 2 + p 22 4P 22 [H 1 "] 4 + p 2222 p 222 p 22 ) 

(23) 

E, = E=[E 3 q .(KH lq .)(K w ) q '[H + ]- q '] 



(24) 



E 2 = EE = [Z 2 q . (KH 2q .) (Kw)'^] ''] 



(25) 

Now, in our case: x = [U0 2 2+ ] and y = [HCit 3- ], therefore 
Eqs. (13-17) lead to: 

C° L = x [Ay + 2A 2 iy 2 ] + 2x 2 Cy 2 + Dy (26) 

C° M = x [1 + Ay + A^y 2 + E] + 2x 2 [EE + Cy 2 ] (27) 

AC° = (C° L - C° M ) = Dy - x[l - A^y 2 + E] - 2x 2 EE (28) 

The variables x and y which are the solutions of Eqs. (26 - 
27), depend on the pH of the medium via the expressions of 
the parameters A, A 2 i, C, D and E, EE. On the other hand, for 
a given total concentrations C° M and C° L and for a given pH, 
the knowledge of x and y allows us to calculate the fraction 
in % of the different species in solution and, therefore, to 
obtain the corresponding diagram of speciation: 

%[M p L n ] = 100.(C°m )-'Pp„ W p [y ] n (29) 

%[M p L n H m ] = 100.(C° M X'Ppnm P P „ [x] p [y ] n (30) 

[H + ] m % [M p L n H m (OH) q ] = 100.(C° M ) _1 P P „ mq P P „ m P P „ 

W p [y] n [H + ] m [H + ]- q (31) 

%[M p (OH) q >] = 100.(C° M r 1 (KH pq .) [OH-] q '[x] p (32) 

V 5 . RESOLUTION OF THE EQUATIONS OF CONSERVATION IN 
THE CASE OF THE URANIUM COMPLEXATION WITH CITRIC 
ACID 

A. Principle 

For a given total concentrations C° M and C° L of M and L, 
the resolution of the system of the conservation equations 
given above depends on five hydrolysis constants KHi q ' , 
three acidity constants K a 1 and on ten complexation constants 
P (see Table 4 and Table 5). According to the literature, these 



last constants are regrouped in seven "sets" which are 
summed up in Table 4. Indeed, this table lists different P 
constants given in the literature [2,3,5,10], or calculated 
according to their corresponding expressions given in "Table 
3". In fact, the examination of this table shows that in reality 
we have only two different independent sets: "Set III" [2,3,5] 
and "Set V" [2,10]. Indeed, the five other sets are some 
particular cases of these two sets obtained by neglecting some 
[MpL n H m (OH) q ] entities (i.e. by taking their corresponding 

Ppnmq 0). 

Note also that, in order to calculate the p 122 constant given 
in "Set V", we used the K 22I constant given in reference [10] 
corresponding to the equilibrium: (U0 2 ) 2+ + 2[H(HCit)] 2 " o 
[(U0 2 )H 2 (HCit) 2 ] 2 - (See Tables 2 and 3). This K 22I constant 
was in fact rejected by the authors of the cited review [10] for 
the following reasons: The original reference doesn't mention 
explicitly the experimental conditions used for the 
determination of this K 22I constant. Moreover, the complex 
[(U0 2 )H 2 (HCit) 2 ] 2 " has not been reported by any other 
researcher. 

On the other hand, Table 5 gives the first five hydrolysis 
constants KH lq . [10] and the three acidity constants K a i 
defined in paragraph 2.1 [7, 10]. All these constants are 
common to the different sets. 

The resolution of the above system involving Eqs. (20 - 32) 
is generally achieved by the numerical methods explained in 
paragraph 3.2. Details are given in Appendix I. However, this 
system could be resolved analytically if we neglect the 
formation of the following species: [(U0 2 )(HCit) 2 ] 4 ~; 
[(U0 2 )H 2 (HCit) 2 ] 2 -; [(UO z ) 2 (OH)] 3+ and [(U0 2 ) 2 (OH) 2 ] 2+ and 
therefore their corresponding constants: (Pi 2 ), (Pi 22 ), KH 2 | and 
KH 22 . In this approximation, both analytical and numerical 
calculations lead to identical speciation diagrams. The main 
advantage of the analytical method consists in the fact that we 
can resolve the system of Eqs. (20-28) by a Microsoft Office 
Excel program in order to calculate x and y for different pH 
and therefore to obtain automatically the graphs of the 
corresponding speciation diagram: i.e. the % [species] in 
function of the pH. Details of the analytical resolution are 
given in Appendix II. It is important to mention that both 
analytical and numerical calculations must verify the 
following mass conservation conditions for (U0 2 ) and (HCit): 

Ipln Em S q p( % [M p L n H m (OH) q ] ) = 100 (31) 

Ipln SmE q n(%[M p L n H m (OH) q ]) = 100 (32) 

However, analytical and numerical calculations show that, 
for some ranges of pH, and even if we use a double numerical 
precision, these conditions are not verified perfectly. The 
reason for that is the very high sensitivity of the calculations 
to the parameters used in these pH ranges. This explains why 
the speciation diagrams presented below are given for 
different domains of pH between and 14. Indeed, for each 
figure, we have chosen the best range of pH for which the Eq. 
(26) and Eq. (27) are verified better than 1%. 

This high sensitivity of the calculations to some 
parameters incited the authors to use their own analytical 
Excel program and their numerical Basic program, in order to 
allow a better control of the different stages of calculations. 
Indeed, this kind of control is not always obvious with the 
utilization of some usual specific software. 
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Note finally that ionic strength corrections can be ignored 
in a first approximation. Indeed more precise numerical 
calculations on the basis of the self consistent procedure led to 
the following conclusion: considering the uncertainties on the 
different constants used in "Set III", "Set IV" and "Set V", the 
improvement of the results by the recurrent method is not 
always significant. Besides, in our calculations and 
experiments, the ionic strength I never exceeded 0.02M. 

TABLE 4 DIFFERENT "SETS" OF COMPLEXATION CONSTANTS OF URANIUM IN 
AUEOUS SOLUTION AT 25°C 



species[(U0 2 ) 2 (OH)] 3+ ,[(U0 2 ) 2 (OH) 2 ] 2+ , [(U0 2 ) (H Cit) 2 ] 4 " 
and [(U0 2 )H 2 (HCit) 2 ] 2 - (i.e. KH 21 = KH 22 = p 12 = Pi 22 = 0). 





"Set 


"Set 


"Set 


"Set 


"Set 


"Set 


"Set 




la" 


lb" 


I" 


II" 


III" 


IV" 


V" 


Pll 


1.0 


1.0 


1.0 


6.3 


1.0 


6.3 


6.3 10 8 




10 7 


10 7 


10 7 


10 8 


10 7 


10 8 




Pill 





7.94 
10 4 


7.94 
10 4 


21.37 


7.94 
10 4 


21.37 


21.37 


pill, 





3.16 


3.16 


4.07 


3.16 


4.07 


4.07 






10 4 


10 4 


10 2 


10 4 


10 2 


10 2 


Pm 








1.26 


1.99 


1.26 


1.99 


1.99 








10" 


10 4 


10" 


10 4 


10 4 


Pl2 














1.0 
10" 


1.0 

10" 


1.0 

10" 


P,22 




















1.0 

10 12 ? 


P22 


1.0 


1.0 


1.0 


7.94 


1.0 


7.94 


7.94 




10 18 


10 18 


10 18 


10 18 


10 18 


10 18 


10 18 


P222 


1 


1 


1 


1 


1 


1 


1 


P2222 





6.3 
10 2 


6.3 
10 2 


79.43 


6.3 
10 2 


79.43 


79.43 


P224 








3.16 


3.89 


3.16 


3.89 


3.89 








10" 


10 10 


10" 


10 10 


10 10 



TABLE 5 HYDROLYSIS CONSTANTS KHlQ' AND ACIDrTY CONSTANTS KA L AT 

25 °C 
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KH21 


KH22 


3.98 10 s 


1.99 10 16 


5.01 10 21 


1.99 10" 


2.51 10 22 


K a , 


K% 


K a _, 






7.94 10 5 


1.99 10 1U 


1.99 10" 







B. Results in Particular Cases 

For a given concentration ratio, Rc = C°m / C°l, the 
analytical method allows us to study easily the sensitivity of 
the diagram to the different constants. We used three 
concentration ratios: Rc= 0.025, Rc= 1 and Rc= 2. Results are 
summarized in "Figs.l to 11". Each figure corresponds to a 
given "set" of constants and to a given ratio R<.. 

1) Results Relative to "Set I" Constants in Aqueous 
Solution at 25°C 

We have first used both analytical and numerical methods 
to calculate the uranyl speciation diagram according to the 
"Set I a " data given in Table 4 (see Fig. 1). This calculation 
corresponds to the simplest case for which we neglected both 
protonation and hydrolysis of M p L n entities. Indeed, "Set Ia" 
differs from "Set I" by the fact that, in this case, all P pnm and 
Ppnmq are ignored. We have also neglected the formation of 





PH 



Fig. 1 Speciation diagram of TJ0 2 (4.10TVI)) in aqueous HCit (2.10 M) at 
25°C with "Set Ia" constants 





PH 



Fig. 2 Speciation diagram of UO 2 2+ (5.10" 4 M) in aqueous HCit 3 (2.10~ 2 M) at 
25°C, with "Set lb" constants and (01 12 = P224 = 0) 




Fig. 3 Speciation diagram of U0 2 (5. 10TVI)) in aqueous HCir (2.10 M) at 
25°C, with "Set I" constants and fJ31 12 > 0, (5224 > 0) 




Fig. 4 Speciation diagram of UO2 (5. 10TVI) in aqueous HCit (2. 10 "M) at 
25°C, with "Set II" constants and (pi 12 > 0, P224 > 0) 
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1 2 3 4 5 6 7 8 9 10 11 12 13 14 



pH 

Fig. 5 Speciation diagram of U0 2 2+ (5.10 3 M) in aqueous HCit 3 " (5.10""M) at 
25°C, with "Set II" constants and (J31 12 > 0, P224 > 0) 

In the second step, we replaced "Set I a " by "Set I b " in 
order to take into account the protonation and the hydrolysis 
of the complexes [(U0 2 )(HCit)]~ and [(U0 2 ) 2 (HCit) 2 ] 2 -; we 
thus obtained the diagram reported in Fig. 2. Comparison 
between Fig. 1 and Fig. 2 shows that, in the pH range: 3 < pH 
< 9, the two diagrams are similar. In both figures we observe 
two plateaus of which the levels are around 35% and 25%. 
The main difference consists in the fact that, in the case of 
Fig. 1 with "Set I a ", they correspond respectively to 
[(U0 2 ) 2 (HCit) 2 ] 2 " (p 22 ) and [(U0 2 )(HCit)]" (p„) species, 
whereas in Fig. 2 with "Set I b ", they correspond respectively 
to the protonated and to the hydrolysed species: 
[(U0 2 ) 2 H 2 (HCit) 2 (OH) 2 ] 2 - (P 2222 ) and [(U0 2 )H(HCit)(OH)r 
(Pi 111)- Besides, around pH «2, we observe in (Fig. 2) the 
apparition of a peak corresponding to a maximum of 
protonation of type [(U0 2 )H(HCit)] (p m ). 

In the third step, we replaced "Set I b " by "Set I" in order to 
take into account the presence of the two species [MLH 2 ] and 
[M 2 L 2 H4] (Pi 12 > 0; p 22 4 > 0). The corresponding diagram for: 
C° M = 5.0 10^M and C° L = 2.0 10" 2 M (!*<;« 1), is given in 
Fig. 3. The comparison between Fig. 2 and Fig. 3 shows that 
the two above-cited species don't modify the shape of the 
diagram qualitatively and quantitatively. The small difference 
consists in the apparition (in the region: < pH < 3) of two 
bell-shaped peaks, with a respectively heights of about 15% 
and 1.5% (at pH a 1.5) corresponding to [MLH 2 ] and 
[M 2 L 2 H4] species. Except this small difference, the two 
figures show: 

The % [U0 2 + ] decreases from 100% at (pH a 0) to 0% at 
(pH«2.5). 

The % [MLH] varies as a bell-shaped curve (in the region: 
< pH < 5); the height of its maximum is of about 85 % at 
(pH a 2.5). Note that this peak is (in our case) non observable 
by capillary electrophoresis method because [U0 2 H 2 Cit] is a 
neutral species. 

In the pH range: 4 < pH < 8.5, both %[MLH(OH)l 
and %[M 2 L 2 H 2 (OH) 2 ] are constants (two plateaus) 
respectively equal to 25% and 35%. At pH a 3.25 and pH a 9, 
the percentage of these two species are sensibly 
equal: % [MLH(OH)] = % [M 2 L 2 H 2 (OH) 2 ] a 20%. 

The %[M(OH) 3 ] presents a gap of height of about 95%, in 
the region: 8.5 < pH < 9.5 

The %[M(OH) 2 ] is a bell-shaped curve (in the pH range: 
8.5 < pH < 10.5); the height of its corresponding peak is about 
13%atpH = 9.3. 
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In summary, according to "Set I" constants and for R; « 
1, we observe quantitatively only 8 species in the pH range: 
< pH < 14. Notice however that, only the two anionic 
complexes: [(U0 2 ) H(HCit)(OH)] _ and [(U0 2 ) 2 H 2 (HCit) 2 
(OH) 2 ] 2 ", are observable by the capillary electrophoresis 
method in the pH range under consideration. Within the limits 
of experimental precision, these two species are observable 
for 3 < pH < 9. 

2) Results Relatively to "Set II" Constants in Aqueous 
Solution at 25°C 

The difference between "Set I" and "Set II" consists 
essentially in the fact that p n (Set I) < Pn(Set II) and p m (Set 
II) « Pin(Set I). The corresponding diagram, for C°m = 
5.0.10" 4 ; C° L = 2.0 10" 2 (Rc« 1), is given in Fig. 4 which 
shows that: 

The peak of the bell-shaped curve (in the region < pH < 
5 of Fig. 3), corresponding to %[MLH], becomes negligible. 
Indeed, its height is only of about 3 % at pH a 2.5. 

A third "plateau" appears, in pH range: 2.5 < pH < 9, with 
a height of about 40 %, which corresponds to the [ML] specie 
(it is the main difference between "Set I" and "Set II"). 

In the same pH range: 2.5 < pH < 9, and by comparison to 
Fig. 3, we observe an inversion of the plateaus corresponding 
to the species [MLH(OH)] (35%) and to the species 
[M 2 L 2 H 2 (OH) 2 ] (10%) rather than the corresponding 
percentages: 25% and 35% respectively observed in the case 
of Fig. 3 relatively to "Set I". 

Except these differences, Fig. 4 is comparable to Fig. 3. 
Notice, however, that in the present case (Fig. 4), the 
intensities of the peaks corresponding to the species [MLH 2 ] 
(30%) and [M 2 L 2 H,] (8%) are the double of those of Fig. 2 
relatively to "Set Ib" 

In summary, with "Set II" constants and Rc« 1, we again 
observe quantitatively 8 species in the pH range: < pH < 14. 
However, in contrast with "set I", three anionic complexes: 
[(U0 2 )(HCit)]-, [(U0 2 )H(HCit)(OH)]- and [(U0 2 ) 2 H 2 (HCit) 2 
(OH) 2 ] 2 ~, are observable by capillary electrophoresis, for 3 < 
pH < 9, within the limits of experimental precision. 

3) Concentration Effect at 25°C 

In order to study the sensitivity of the speciation diagram 
to the concentration ratio, R; , in the case of the "Set II" 
constants, we established the diagrams presented in Fig. 5 and 
Fig. 6 respectively: 

C° M = 5.0 10" 4 ; C° L = 5.0 10 4 ^ Re = 1, (given in Fig. 5). 

C° M = 1.0 10" 3 ; C° L = 5.0 10" 4 <^ Rc= 2, (given in Fig. 6). 

Fig. 5 shows that, when the solution is equimolar in 
[(U0 2 )] and [(HCit)], peaks corresponding to species [MLH], 
[MLH 2 ] and [M 2 L 2 Ht] are negligible (< 2%). In contrast, the 
height of the [M(OH) 2 ] peak (in the region 7 < pH < 10) is 
now of about 30 % at pH a 9. 

However, and as in the case of Fig. 4 (Rc«l), Fig. 5 
shows three plateaus corresponding to the % of the species: 
[ML], [MLH(OH)] and [M 2 L 2 H 2 (OH) 2 ], of heights 
respectively equal to 40%, 35% and 10%, in the pH range: 3 < 
pH < 8. Notice also that the two gaps corresponding to the 
species [(UO z )] 2+ and [(UO z ) (OH) 3 ] _ (for R, = 1) are 
comparable to those of Fig. 4 (for R; «1). 
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In conclusion, we can say that, for R c =l, the species 
([(U0 2 )(HCit)r, [(U0 2 )H(HCit)(OH)]" and 

[(U0 2 ) 2 H 2 (HCit) 2 (OH) 2 ] 2 ") are observable by capillary 
electrophoresis, in the pH range: 3 < pH < 8 . 




0,5 1 1,5 2 2,5 3 3,5 4 4,5 
pH 



Fig. 6 Speciation diagram of U0 2 2+ (10 3 M) in aqueous HCit 3 " (5.10" 4 M) at 
25°C, with "Set II" constants and ((31 12 > 0, 0224 > 0) 
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Fig. 7 Speciation diagram of U0 2 2+ (4.10 4 M)) in aqueous HCit 3 (2.10 2 M) at 
25°C with "Set III" constants; pll2 = 




Fig. 8 Speciation diagram of U0 2 2+ (4.10^)) in aqueous HCit 3 (2.10 2 M) at 
25°C with "Set IV" constants; pi 12 = 

On the other hand, Fig. 6 (for Rc = 2) shows that the 
excess of [(U0 2 )] 2+ ions changes appreciably the speciation 
diagram in the sense that the respective heights of the plateaus: 
[ML], [MLH(OH)] and [M 2 L 2 H 2 (OH) 2 ] are decreased, in the 
pH range: 3 < pH < 5, toward: 15%, 10% and 5%, 
respectively. Notice also that the variation with pH of the 
percentage of the species [(U0 2 )] 2+ is now completely 
different, by comparison to Figs. 4, 5. In particular, this 
percentage remains sensibly constant equal to 45 % in the pH 
range: 3 < pH < 5 (for Rc = 1). Nevertheless, we can say that, 
for R-2, the species [(U0 2 )(HCit)] _ and [(UO z ) H(HCit) 
(OH)]" remain observable by capillary electrophoresis, in the 
pH range: 3 < pH < 5. However, and within the limits of the 
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experimental precision, this fact is less true for 
[(U0 2 ) 2 H 2 (HCit) 2 (OH) 2 ] 2 - species. 

C. Results in the General Case at 25°C 

In the general case, and in addition to the previous entities, 
we have taken into account the formation of the following 
species: [(U0 2 )(HCit) 2 ] 4 ": (p 12 ), [(U0 2 ) 2 (OH)] 3+ : (KH 21 ), 
[(U0 2 ) 2 (OH) 2 ] 2+ : (KH 22 ). Their corresponding constants: (p 12 ), 
(KH 21 ) and (KH 22 ) are given in Table 4 ("Set III" = "Set I" + 
P12 or "Set IV" = "Set II" + p 12 ) and in Table 5. The 
questionable case of the specie [(U0 2 )H 2 (HCit) 2 ] 2 - : (p l22 ) will 
be analyzed at the end of this paragraph. 

The numerical result corresponding to "Set III" is 
represented in Fig. 7. Comparison with Fig. 3 ("Set I") shows 
that except for the two pH ranges: 3 < pH < 5 and 9 < pH < 11, 
the two diagrams are completely different. Indeed, with "Set 
III" constants and for pH > 5, the percentages of both species: 
[(U0 2 )H(HCit)(OH)]": (p lm ) and [(U0 2 ) 2 H 2 (HCit) 2 (OH) 2 f: 
(P2222 ) decrease with pH from 30% to respectively: 10% and 
5% in the pH range: 6 < pH < 9. This decrease is offset by the 
apparition in the same pH range of an important "plateau" of 
about 80 % height, corresponding to the specie 
[(U0 2 )(HCit) 2 ] 4 -: (p 12 ). 

Now if we replace "Set III" by "Set IV", so that p u (Set III) 
< Pn(Set IV) and p m (Set IV) « p m (Set III) (by analogy to 
the transformation of "Set I" into "Set II"), we obtain the 
diagram given by Fig. 8. Comparison of Fig. 8 with the 
previously Fig. 7 indicates the apparition with a relatively 
important percentage of the specie [(U0 2 )(HCit)]" in addition 
to the species: [(U0 2 )H(HCit)(OH)]" and [(U0 2 ) 2 H 2 (HCit) 2 
(OH) 2 ] 2 " (by analogy to the comparison between Fig. 4 ("Set 
II") and Fig. 3 ("Set I")). Consequently, the height of the 
"plateau" corresponding to the specie [(U0 2 )(HCit) 2 ] 4 ~is now 
of about 60 % rather than 80%. 

In conclusion, we can say that numerical calculation for 
Re«l predicts that the species [(U0 2 )(HCit)]~, 
[(U0 2 )H(HCit)(OH)]- and [(U0 2 ) 2 H 2 (HCit) 2 (OH) 2 ] 2 - should be 
observable by capillary electrophoresis in the region: 3 < pH < 
5, even if we take into account of the species [(U0 2 )(HCit) 2 ] 4 ~ 
which become observable only for pH>5 and this 
independently of the choice of the literature values of Pnand 
Pin for the calculation. Indeed, Pnand p m influence rather 
the relative peak intensities of the three observable entities. 

If in addition we suppose the formation of the specie 
[(U0 2 )H 2 (HCit) 2 ] 2 in conformity with the "Set V" = ["Set IV" 
with: Pi 22 =1.0 10 12 ] we therefore obtain the diagram given in 
Fig. 9 calculated for Rc«l. This diagram shows clearly that 
the introduction of the Pi 22 constant provokes the complete 
disappearance of all the species [(U0 2 )(HCit)]", [(U0 2 ) 
H(HCit)(OH)]and [(U0 2 )2H 2 (HCit) 2 (OH) 2 ] 2 -in the pH range: 
2 < pH < 5. This disappearance is compensated (up to 100% !) 
by the species[(U0 2 )H 2 (HCit) 2 ] 2 " . For pH>5, Fig. 9 remains 
quite similar to Fig. 8. Thus, according to the above 
hypothesis (p 12 2 *0), numerical calculations show that for 
Rc«l, [(U0 2 )H 2 (HCit) 2 ] 2 ~ is the only observable specie in the 
pH range: 3 < pH < 5. We will see that this conclusion is not 
verified by the capillary electrophoresis experiments. 

1 ) Concentration Effect in the General Case at 25°C 

In order to study, in the general case characterized by 
"Set V" constants, the sensitivity of the speciation diagram to 
the concentration ratio Rc, we have established the following 
diagrams given by Fig 10 and Fig. 1 1, respectively for: 
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C°m=5.0 10- 4 ;C° L = 5.0 10 _4 <^R C = 1. 

C° M = 1.0 10" 3 ; C° L = 5.0 10" 4 <^ Rc= 2. 

Comparison of Fig 10 (with: p 12 > 0, Pi 2 2> 0) with Fig. 5 
(P12 = 0, P122 = 0) shows that the formation of species 
[(U0 2 )(HCit) 2 ] 4 " has no significant influence on the diagram, 
and that the influence of [(U0 2 )H 2 (HCit)2] 2 " species is also 
negligible by comparison with Fig. 9 (Rc « 1). Indeed, the 
height of their corresponding bell-shaped peak in Fig. 10 
doesn't exceed 7% in the pH range: 2 < pH < 4. 
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Fig. 9 Speciation diagram of U0 2 2+ (4.10 4 M)) in aqueous HCit 3 (2.10 2 M) at 
25°C with "Set V" constants; pi 12 > 




1Q -0-%(U02)2(OH) 
^%(U0212(OH)2 



Fig. 1 Speciation diagram of U0 2 2+ (5 . 1 O^M) in aqueous HCit 3 (5 . 1 4 M) at 
25°C with "Set V" constants; (31 12 > 
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Fig. 1 1 Speciation diagram of U0 2 2+ (10" 3 M) in aqueous HCit 3 (5.10" 4 M) at 
25°C with "Set V" constants; (31 12 > 



Note also 

n 2+ 



that the two gaps corresponding to both 
[(U0 2 )] z+ and [(U0 2 )(OH) 3 ]" species, for Rc = 1, are also 
comparable to those of Fig. 5 (Rc =1). 



In the same manner, comparison of Fig. 1 1 (R c = 2, Pi 2 > 0, 
P122 > 0) with Fig. 6 (relative to "Set II": p 12 = 0, Pi 22 = and 
Rc = 2), shows that the formation of the two species 
[(U0 2 )H 2 (HCit) 2 ] 2 ~ and [(U0 2 )(HCit) 2 ] 4 ~ is always negligible. 
Indeed, the heights of their corresponding peaks are 



respectively less or of about 3% in the pH range: 2 < pH < 3. 
Notice finally that, according to Fig. 11, the relative excess of 
(U0 2 ) 2+ ions (characterized by the ratio: Rc = 2), increases the 
maxima of the bell-shaped peaks relative to the species 
[(U0 2 ) 2 (OH) 2 ] 2+ and [(U0 2 )(OH) 2 ], respectively up to 43% 
in the pH range: 3 < pH < 8, and to 40% in the pH range: 6 < 
pH < 10. However, these peaks cannot be observed 
experimentally in the present case of our experience. 

In conclusion, the formation of [(U0 2 )H 2 (HCit)2] 2 ": (p 122 ) 
and [(U02)(HCit) 2 ] 4 ": (P12) species seems to appear with 
important and constant percentages (respectively: 95% for: 
1.5 < pH < 6.5 and 60% for: 4.5 < pH < 10), only for low 
ratios: Rc « 1. In contrast, for Rc > 1, the percentage of 
[(U0 2 )H 2 (HCit) 2 ] 2 - is < 5%, whereas, the rate of [(U0 2 ) 
(HCit) 2 ] 4 " remains negligible. This means that, in the range of 
pH studied experimentally (pH < 5), numerical computations 
predict the apparition of only one intense peak for dilute 
solutions in uranium which must disappear with increasing 
the U (VI) concentration. This prediction is not in agreement 
with our previous observations [16]. Therefore, the hypothesis 
of the existence of the [(U0 2 )H 2 (HCit) 2 ] 2 " specie seems 
unlikely. In contrast, calculations according to: "Set I to Set 
IV" predict the observation of three peaks in the pH range: 2 
< pH < 5, corresponding to the [(U0 2 )(HCit)]\ [(U0 2 ) 
H(HCit)(OH)]- and [(U0 2 ) 2 H 2 (HCit) 2 (OH)2] 2 " species. 
However, the importance of the [(U0 2 )(HCit)]" peak depends 
on the choice of the p n and p m constants. Notice that 
according to the reference [10], the most probable species in 
the pH range: 2 < pH < 5 are rather [(U0 2 )(HCit)]" and 
[(U0 2 ) 2 (HCit) 2 r (See Fig. 1 obtained with "Set la"). This 
conclusion is not compatible with the diagrams calculated 
with the data of the all other sets as long as we assume that 
the different protonation and hydrolysis reactions have 
reached their complete chemical equilibrium. 
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Fig. 12 Comparison between electropherograms of sodium citrate buffer 
solutions at pH=3 containing or not NO3" ions 
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Fig. 13a Electropherograms of 1.1(T 5 M nitrate and 5.10 4 M U(VI) in 0.02 M 
citrate buffer at pH=2.5 and pH=3.0 
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Fig. 13b Electropherograms of 1.10 5 M nitrate and 5.10" 4 M U(VI) in 0,02M 
citrate buffer at pH=3.5 and pH=4.0 
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Fig. 14 Electropherograms of 5.10 4 M U(VI) in 0.02M citrate buffer at 
various pH values 

VI EXPERIMENTAL 

A. Apparatus 

The capillary electrophoresis apparatus is a modular 
system consisting of a Spectraphoresis 100 injector 
(hydrodynamic mode) coupled with a high voltage (0-30kV) 
power supply (Prime Vision VIII from Europhor) and a 
scanning UV-visible detector (Prime Vision IV from 
Europhor). In the present study, a voltage of 30kV was 
applied and a fixed wavelength of 200 nm was chosen for the 
absorbance measurements. The external temperature was 
maintained equal to 25°C. The acquisition of the 
electropherograms and the data treatment were performed 
using the chromatography software BORWIN (JMBS 
Developpements). 

Each measurement was repeated two or more times in 
order to assure the reproducibility condition. 

B. Chemicals 

All solutions were prepared with deionised water and 
chemicals of analytical reagent grade. For the sodium citrate 
buffer solutions, we have used different Fluka solutions for 
HPCE (210" 2 M) having pH values between 2.5 and 5.5. The 
uranyl stock solution (concentrated U0 2 (C10 4 ) 2 ) was 
prepared by dissolving lg of uranyl nitrate hexahydrate salt, 
U0 2 (N0 3 )2.6H 2 (FLUKA), in 25 ml of concentrated 
perchloric acid. The resulting mixture was then almost 
completely evaporated at 150°C in a sand bath and this 
procedure was repeated several times in order to eliminate all 
the nitrates. This acidic (HC10 4 + U0 2 (C10 4 ) 2 ) stock solution 
was then diluted in order to adjust the uranyl concentration to 
around 10" 2 M (the exact concentration was determined by 
fluorescence spectroscopy). 
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For each pH value of the citrate buffer, a corresponding 
sample was prepared by mixing appropriate volumes of 
diluted uranyl stock solution and sodium citrate buffer 
solution. The final pH of this sample solution was adjusted by 
adding a small quantity of NaOH solution (0.1M). The total 
final volume of the different prepared samples is equal to 
1000 |al. 

It is important to mention that, before performing any 
electrophoresis experiment, the sample solutions were kept at 
rest for at least two days, in order to allow the different 
reactions of complexation, protonation, and hydrolysis, to 
reach their complete chemical equilibrium. 

C. Procedure 

The silica capillary (SUPELCO, 75 |im of internal 
diameter, 75 cm long), was first conditioned by successive 
washes with 0.1 NaOH, then by deionised water, and finally 
by 0.02M sodium citrate buffer solution of a given pH. 

The cleanliness of the capillary is checked by a first run 
using only the buffer solution under study (check of the 
baseline). The corresponding electropherogram should not 
present any peak but only a smooth baseline (or noise line). 

In order to control the eventual effect of the 
electroosmotic flow (neglected at low pH values [11]), as well 
as other experimental fluctuations, which can shift the peak 
positions of the anionic species under observation, nitrate ions 
(at very low concentrations) were added to some samples. 
This anion has already been used in previous works as a 
reference anion for mobility measurements [12, 13, 16]. Fig. 
12 gives an example of such electropherograms obtained at 
pH = 3 (N0 3 " peak appears around 2 min). 

Each run consisted in a hydrodynamic injection of 2n 
nanoliters during n sec, of a sample solution (uranyl solution 
in buffer citrate with or without nitrate ions) into the capillary. 
The corresponding electropherogram was obtained by 
applying 30 kV during an acquisition time of at least 10 min. 
Finally, the capillary was carefully rinsed with the buffer 
under study between two runs and kept filled with deionised 
water during the night. 

D. Results and Discussion 

As expected, since citrate ions are absorbing species in the 
UV domain, the stability of the baseline is very sensitive to 
the buffer pH value, as can be seen in Figs. 13, 14. The 
background absorbance was observed to increase with pH and 
strong fluctuations to appear for pH > 4.5. These changes in 
the level and quality of the baseline can be assigned to an 
augmentation of the charged citrate ions of the buffer when 
pH increases. Consequently, the current passing through the 
capillary increases sharply with pH: (from 20 □ A at pH 2.5~ 
3 to 135 U A at pH 5). For these reasons, the best results were 
obtained, for: 2.5 < pH <4.5 (the reasonable limit, 70 - 75 DA, 
is reached at pH 4.5). The different results are summarized in 
the electropherograms presented in "Figs. 13a, 13b" and in Fig. 
14. All of them were obtained for C° L = 0.02 M, C° M = 5 10" 
4 M and eventually C° N0 3- = 1.10" 5 M. They show the existence 
of three peaks which can be assigned with the anion uranyl 
species; the most mobile entities are those which present, the 
low time of migration. 

1 ) Electropherograms for: 2 < pH <3 
According to "Set I" constants (p n , p lu , p lm , (3 U2 , p 22 , 
P2222, P224, KH n , KH 12 , KH 13 ) or "Set III" = {"Set I" + p 12 ) 



IJNESE Vol. 2 No. 1 2012 PP. 11-22 www.ijnese.org ©World Academic Publishing 

- 19- 



The International Journal of Nuclear Energy Science and Engineering 



IJNESE 



(see "Fig.3" or "Fig. 7"), only cationic ([(UO z )H 2 (HCit)] + , 
[(U0 2 ) 2 H4(HCit) 2 ] 2+ ) or neutral ([(U0 2 )H(HCit)]) species, 
exist for pH < 3. The two anionic species: ([(U0 2 ) 
H(HCit)(OH)]-, [(U0 2 ) 2 H 2 (HCit) 2 (OH) 2 ] 2 ), are thus expected 
to appear quantitatively only for pH >3.5. Notice that, 
according to "Set II" constants (P'n, P'm, P'mi, P'm, P'22, 
P'2222, P'224, KH„, KH 12 , KH„) or "Set IV" = {"Set II" + p 12 } 
constants (see Fig.4 or Fig. 8), the three anionic species: 
([(U0 2 )(HCit)r, [(U0 2 )H(HCit)(OH)]- and [(U0 2 ) 2 H 2 (HCit) 2 
(OH) 2 ] 2 ~) appear also quantitatively, only for pH > 3. These 
theoretical results are in good agreement with the fact that, for 
pH = 2.5 and pH = 3, the corresponding electropherograms do 
not point out any peak associated with anionic uranyl species. 
The only apparent peak shown in Fig. 13a, with a migration 
time t « 2 min, is due to the added nitrate ions. 

2) Electropherograms for: 3 < pH < 5 

In contrast, for buffer pH values around 4, the 
corresponding electropherograms show three separate peaks 
at migration times of * 3.2 min, t « 4.7 min and t « 5 min (see 
Fig. 13b, 14). Owing to our theoretical results, these peaks 
could be assigned to the three following species: 
[(U0 2 ) 2 H 2 (HCit) 2 (OH) 2 ] 2 ", [(U0 2 )H(HCit)(OH)]" and [(U0 2 ) 

(Hat)]-. 

Their electrophoretic mobility, u ; , can be derived from the 
one of NCV: Umm, by applying the following expression: 



Ui = u N(M + (d ITU) (1/tMi - l/t N o3 _ ) 



(33) 



where u N0 3 is the ionic mobility of NCV (7,4 10~ 4 cm 2 V" 1 
s" 1 ); U is the applied voltage (30 000 V); L is the capillary 
length (75cm) and d, the length from the capillary inlet to the 
detector (35 cm). The parameters: t M i and t N0 3 are respectively 
the times taken by the ionic "i" specie and by the NO3" specie 
to reach the detector [12]. 

TABLE 6 ELECTROPHORETIC MOBILITIES, UI, CORRESPINDING TO THE PEAKS 
OBSERVED IN FIG. 13B (N03- BEING TAKEN AS REFERENCE) 





pH = 3.5 




pH = 4.0 






t(migration) 

(sec) 


10 4 Ui 
(cm 2 V" 1 


t(migration) 
(sec) 


10 4 Ui 
(cm 2 V" 1 
s" 1 ) 


no 3 " 


122.2 


(7.4) 


128 


(7.4) 


species 
1 


180.9 


5.08 


195 


5.06 


specie 2 


256 


3.66 


279.5 


3.73 


specie 3 






294.8 


3.54 



The calculation leads to the following electrophoretic 
mobility values: 5.1 W 4 , 3.7 10" 4 and 3.5 10" 4 cm 2 V" f s 1 (see 
Table 6), corresponding to the three uranyl complex forms 
detected in Fig. 13b. However, the exact correspondence 
between each of these values and the above cited anionic 
species is not a priori evident. Nevertheless, the ionic 
mobilities of monomer species are generally assumed to be 
approximately proportional to their charge over size ratio [13]. 
Indeed, according to the Nernst -Einstein relation: UjDj" = 
(ZieXkT)" 1 and Stokes-Einstein relation: D, = kr^^R) 1 , 
where: D ; is the self diffusion coefficient, kT is the thermal 
energy, r\ is the viscosity and Z,e, R; are the charge and the 
radius of the considered "i" specie, the expression of the 
mobility is therefore: U; = (Ze)/ (67rnR). Consequently, we 
have attempted to make the following qualitative predictions 
on the basis of the different structures attributed by recent 
literature to the three complexes under consideration: 



[(U0 2 )H(HCit)(OH)]- and [(U0 2 )(HCit)]" they have the 
same charge and similar radii; therefore, they should have 
near mobility values. The two peaks observed at pH > 3.5 at 
migration time t « 4.8 min, and characterized by a weak 
separation and a comparable height(see Fig. 13b, 14), can thus 
be attributed to these two species. Indeed, if we assume the 
validity of the "Set II" constants, these two species are also 
expected to have close complexation rates (or %), in 
accordance with the theoretical predictions: 40% and 35%, for 
respectively [(U0 2 )(HCit)]" and [(U0 2 )H(HCit)(OH)]" (see 
Figs. 4, 8). We can note, however, that at pH = 3.5, and in 
contrast with pH =4, only one of the two close peaks appears 
clearly on the corresponding electropherograms (the second 
peak could only be guessed on the electropherogram of Fig. 
14). Comparison between Fig. 13a and Fig. 13b shows that 
the percentage of the two above cited species increases from 
pH «2.5 to pH 4 with a gap at pH « 3. This tendency is 
reflected in the diagrams given in Figs. 4, 8 but with a gap at 
pIU 2.5. 

The [(UO z ) 2 H 2 (HCit) 2 (OH) 2 ] 2 - the complexion has a 
double negative charge (Z, = -2) and it exhibits three possible 
conformations [14]: two of them are rather spherical 
characterized by a large radius R ; and therefore by a weak 
diffusion coefficient: D;; In contrast, the third most probable 
conformation, is a cylindrical chain presenting a weak viscous 
friction and therefore, a large mobility u ; [15]. However, since 
its corresponding complexation rate is about 10%, (according 
to our theoretical calculations: Figs. 4, 8), its relative 
electrophoretic peak should be weak. This fact is confirmed 
by the electropherograms performed at pH 3.5, pH4 and pH 
4.5 after a migration time of about: t « 3.0-3.5 min. 

Note moreover that, for pH > 4.5, this peak (at t « 3 min) 
is always observable, whereas the other two peaks are, more 
or less hidden by the noise of the baseline (see Fig. 14). On the 
other hand, according to our calculation using "Set II" 
constants, the increasing of the U(VT) concentration (i.e. the 
ratio Rc), will not give more information concerning these 
three species since their percentages decrease with increasing 
the ratio Rc. 

It is interesting to underline at this stage, and more 
particularly for the speciation measurements under 
consideration, the interest of the capillary electrophoresis 
method by comparison to some other analytical methods: the 
latter ones are rather sensitive to the coordination state of 
U(VI), (for example, because of the variation of the U-O 
distance with pH). It ensues that these techniques cannot in 
general distinguish between protonated and non- protonated 
complexes: M,Li, MiLiH, , M,LiH 2 , M 2 L 2 , M 2 L 2 H 2 , and 
M 2 L 2 H(, as well as between hydrolyzed and non-hydrolyzed 
complexes: (for example: MiLiH^OH), or M 2 L 2 H 2 (OH) 2 ). In 
other terms, and for the present study, they could not inform 
us sufficiently on the protonation state nor on the hydrolysis 
state of M = U(VT). On the other hand, the capillary 
electrophoresis method, which is very sensitive to both the 
charge Z;e and/or the size of the species, allows the separation 
of the different anionic complexes expected to be formed on 
the base of their relative mobility. 

VII CONCLUSION 

In order to complete more deeply our previous work [16] 
in selecting the more appropriate complexation constants 
among those published by the literature, and following the 
controversy concerning the processes of protonation and the 
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hydrolysis of the uranyl-citrate species, we undertook to 
compare our new uranyl-citrate speciation diagrams 
calculated analytically or numerically with experimental 
results obtained by the capillary electrophoresis method, in 
dilute citrate aqueous solutions (0.02M) at 25°C and for: 2.5 < 
pH<5. 

Calculations show 16 possible [M p L n H m (OH) q ] entities, 
with: M = (U0 2 ) 2+ , L = (HCit) 3 \ For R, = 0.02, cationic 
([(U0 2 )H 2 (HCit)] + , [(U0 2 ) 2 H4(HCit) 2 ] 2+ ) or neutral ([(U0 2 ) 
H(HCit)]) protonated U(VI)-citrate-complexes appear for pH 
< 3. Using "Set II" constants, three anionic ([(U0 2 )(HCit)]~ 
,[(U0 2 )H(HCit)(OH)]-, [(U0 2 ) 2 H 2 (HCit) 2 (OH) 2 f) species 
appear mainly for: 3 < pH < 5 with respectively rates of about: 
40%, 35% and 10%. For: 5 < pH < 9, these proportions 
decrease of half (20%, 15% and 5%) if we assume the 
formation of the [(U0 2 )(HCit) 2 ] 4 ~ specie. For pH > 9, the 
specie [(U0 2 )(OH) 3 ]~ is predominant. Notice that if we 
assume the existence of the [(U0 2 )H 2 (HCit) 2 ] 2 ~ specie 
(reported by only one reference), the speciation diagram is 
then completely modified in the range: 1.5 < pH < 6, since its 
percentage in this domain is evaluated to 95%. However, for 
Re > 1, both [(U0 2 )(HCit) 2 ] 4 - and [(U0 2 )H 2 (HCit) 2 ] 2 - species 
become negligible; consequently, the respective proportions 
of the three [(U0 2 )(HCit)]\ [(U0 2 )H(HCit)(OH)]-, [(UO z ) 2 
H 2 (HCit) 2 (OH) 2 ] 2 " species become constant (3 plateaus) in the 
region: 3 < pH < 8. For: Rc = 1, these percentages are about: 
35%, 30% and 10%. For: R c = 2, they decrease to about: 20%, 
17% and 5% because of the predominance of the free (U0 2 ) 2+ 
ion for pH < 4, and also due to the apparition of the 
hydrolyzed ions: [(U0 2 ) 2 (OH) 2 ] 2+ for 4.5 < pH < 7 and 
[(U0 2 )(OH) 2 ] for 7 < pH < 9. It is interesting to notice that 
less recent constant data ("Set lb"), lead to the same results, 
except that the specie [(U0 2 )(HCit)]" is substituted by the 
neutral specie [(U0 2 )H(HCit)], which appears only for pH < 
3.5. 

The application of the capillary electrophoresis technique 
to the study of the uranyl-citrate speciation in dilute acidic 
aqueous solutions for: (2.5 < pH < 5), confirms qualitatively 
our theoretical predictions. Indeed, the corresponding 
electropherograms, show two close electrophoretic peaks (or 
nearly superimposed) of comparable heights and a third small 
peak. We could assign the two first peaks to [(U0 2 )(HCit)]" 
and [(U0 2 )H(HCit)(OH)]" species, while the third peak 
represents the [(U0 2 ) 2 H 2 (HCit) 2 (OH) 2 ] 2 ~ species. The 
corresponding proportions are about 40%, 35% and 10% at 
pH 4. On the other hand, our experimental results invalidate, 
as expected, the existence of the [(U0 2 )H 2 (HCit) 2 ] 2 ~ specie 
(which, in principle, should appear with a rate of about 95% 
for 1.5 < pH < 6.5 and for R;«l). Moreover, the observation 
of the [(U0 2 )(HCit) 2 ] 4 ~ specie in the region of pH > 4.5, seems 
to be screened by the noise effect. Notice however that, this 
numerical speciation analysis is based on the assumption: that 
all the protonated and hydrolyzed species have reached their 
complete chemical equilibrium. Nevertheless, the present 
study enabled us to emphasize the capacity of the capillary 
electrophoresis technique to separate (at least qualitatively) 
different ionic complexes and to specify their degrees of 
protonation and hydrolysis even out of their equilibrium states. 
Indeed, it was particularly possible for us to prove the 
reliability of the processes of protonation and hydrolysis of 
both the monomer and the dimer: [(U0 2 )(HCit)]~, 
[(U0 2 ) 2 (HCit) 2 ] 2 - in the following forms: [(U0 2 ) H(HCit) 
(OH)]" and [(U0 2 ) 2 H 2 (HCit) 2 (OH) 2 ] 2 \ This conclusion opens 
new fundamental prospects in the aim of explaining (in terms 
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of competition between ligands, OH" and H 2 [17]) how the 
complexation of a metal cation similar to U0 2 2+ (whose 
hydrolysis generally occurs at relatively elevated pH) 
facilitates in certain cases the process of its hydrolysis at a 
lower pH. 
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APPENDIX I NUMERICAL RESOLUTION OF THE EQUATIONS OF 
CONSERVATION 

In the general case of: pi 2 >0 and p [22 >0, the 
corresponding system of conservation equations is given by 
Eqs. (26-28) with: x = [UO z 2+ ] and y = [HCit 3 ], therefore: 

y=0.5[-D+A l/2 ](xA 21 )" 1 (1.1) 

with: A = D 2 + 4(xA 2 ,)[AC + x(l+E) + 2x 2 EE] (1.2) 

AC° = (C° L - C° M ) = Dy - x [1 - A 21 y 2 + E] - 2x 2 EE (1.3) 

C° L = x [Ay + 2A 2 iy 2 ] + 2x 2 Cy 2 + Dy (1.4) 

C° M = x [1 + Ay + A^y 2 + E] + 2x 2 [EE + Cy 2 ] (1.5) 

Recall that: x, y, A, A 2 i, C, D, E, EE are functions of the 

pH. 

For a given pH, the numerical method consists to start 
from an initial x value. x can be chosen equal to the final 
solution x' of the above system (1. 1-1.5) obtained for the 
previous pH' = pH - ApH. After n iterations, we define: 

x (n) = x ± nAx. Thus: y(n) = = 0.5[-D(n) +A(n) 1/2 ](x(n)A 21 )" 

(1.6) 

The convergence is reached when both x(n), y(n) values 
verify Eqs. (1.4 and 1.5); within a given precision: %Px = 
100[x(n) - x(n-l)]/x(n). If it is not the case, we proceed in the 
same manner with the next x(n+l), y(n+l) values. It is 
obvious that the precision %Px of these solutions is increased 
when the step Ax is very small. However, this precision is 
limited by the precision of calculation of the computer. In 
general, our convergence precision is less than 1% which is 
comparable to the estimated error on the concentrations. 

APPENDIX II ANALYTICAL RESOLUTION OF THE EQUATIONS OF 
CONSERVATION 

The system of the conservation equations given in 
paragraph 4.3 could be reduced to an algebraic equation in 
"x" of degree 4, if we ignore the following species: 
[(U0 2 ) 2 (OH)] 3+ , [(U0 2 ) 2 (OH) 2 ] 2+ [(U0 2 )(HCit) 2 ] 4 " (p 12 = 0) 
and [(U0 2 )H 2 (HCit) 2 ] 2 " (p 122 = 0). Therefore, EE = A 21 = 0. 
With such an assumption, the system of the conservation 



equations becomes: 

C° L = Axy + 2x 2 Cy 2 + Dy (II. 1 ) 

C° M = x[l+Ay+E] + 2x 2 Cy 2 (11.2) 

AC° = (C° L - C° M ) = Dy - x[l+ E] (JX3) 

Therefore, y = (AC° + x[l+ E])D _1 (JX4) 

And: x 4 + a,x 3 + a 2 x 2 + a 3 x + a 4 = (H5) 
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a 2 


a 3 




2AC°(1+E)-' 


D" 2 [2C(AC°) 2 +A(1 + 
E)D](ao)-' 


[1 + E+A(AC°)D 
'](ao) J 


(- 

C° M )(ao)"' 



With a = 2C(D 2 )(1+E) 2 

The previous Eq. (H5) can be factorized as following: 
x 4 + ajx 3 + a 2 x 2 + a 3 x + a 4 
= (x 2 + xS m +P)(x 2 + xS' m + P') = (11.6) 



After identification we obtain: 

(S m +S' m )=a 1 ; (S m .S' m + P + P') = a 2 ; (P'.S m + PS' m ) = a 3 ; 

P.P' = a 4 (H.7) 

If now we define the parameter "z" as: 

z = (P + P') (H.8) 

We can therefore express the parameters S m , S' m , P and P' 
in terms of "z" and the "a ; , a 2 , a 4 " coefficients: 

S m = 0.5[z-(z 2 -4a 4 ) 1/2 ] ; S' m = 0.5[z + (z 2 - 4a 4 ) 1/2 ] (H.9) 

P = 0.5[a 1+ {(a 1 ) 2 -4a 2 +4z} 1/2 ] ; 

P' = 0.5[a r {( ai ) 2 - 4a 2 + 4z} 1/2 ] (11.10) 

Now, the condition on the "a 3 " coefficient (P'.S m + PS'm) 
= a 3 shows that "z" is solution of the following equation (of 
third degree) in z: 

z 3 + A ; z 2 + A 2 z + A 3 =0 (11.11) 

With: Ai= - a 2 ; A 2 = (a ; a 3 - 4a 4 ) ; 

A 3 =[4a 2a4 -(a 3 ) 2 - ( ai ) 2 a 4 ] (11.12) 

The "z" real solutions depend on the sign of the 
determinant D t : 

D t = Q 3 + R 2 ;Q=[3A 2 -(A 1 ) 2 ]/9; 

R = [9A[A 2 -27 A 3 -2(A!) 3 ]/54 (11.13) 



D, >0 


D, = 


D, <0 


z=S+T- 
A,/3 


z = (S + T)/2 


z = 2(-Q)" 2 Cos{(6/3) + n;t/3)-A 1 /3 ; 
(n = 0, 2, 4) 
Cos(9) = ±R.(-Q 3 )" 2 



with: S=[R+(D) 1/2 ] 1/3 and T= [R - (D) 1/2 ] 1/3 (11.14) 



Finally, the four "x" solutions relatively to the Eq. (H5) 

are: 
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x 1+ = 0.5[-S m + (A)" 2 ] ; Xl . = 0.5[-S m -(A)" 2 ] (11.15) 

x' 1+ =0.5[-S' m + (A') 1/2 ] ;x',. = 0.5[-S' m -(A') 1/2 ] (11.16) 

with: A = [(S m ) 2 - 4P] 1/2 : A' = [(S' m ) 2 - 4P'] 1/2 (11.17) 

Only one solution "x" is physically correct, because x = 
[U0 2 2+ ] = [M] must be positive and must verify the condition 
of conservation of the total mass of both U0 2 and HCit. The 
corresponding "y" value is: 

y = [HCit 3- ] = [L] = (AC° M + x[l+ E])/D (11.18) 
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